Optical Gap Solitons in Nonresonant Quadratic Media 
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We demonstrate an important role of the process of optical rectification in the theory of nonlinear 
wave propagation in quadratically nonlinear [or x ] periodic optical media. We derive a novel 
physical model for gap solitons in x' 2 ^ nonlinear Bragg gratings. 
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As has been recently demonstrated, large optical non- 
linearities can be generated in noncentrosymmetric me- 
dia by means of the so-called cascading effects, due to 
wave mixing under the condition of nearly phase-matched 
second-harmonic generation and other parametric pro- 
cesses §. It has been also shown that cascaded nonlin- 
earities can support spatial optical solitons || and also 
different types of gap solitons in periodic Bragg gratings 
with a quadratic [or x'* 2 ''] nonlinear response 

Importantly, when an input electromagnetic wave E at 
frequency to is launched into a noncentrosymmetric ma- 
terial, it generates also a quasi-static electric field (or dc 
field) at frequency zero. This effect is known as optical 
rectification, and it is usually described by a contribu- 
tion to the medium nonlinear polarization P of the form 
p i = eoXijk(0;cJ, -uj)Ej(ui)E^(ui), where Xijk(0;u,-u) 
is the nonlinear optical susceptibility describing optical 
rectification p|. Such an induced dc field changes a re- 
fractive index via the linear electro-optic effect. As has 
been recently shown by Bosshard et al. || , both theoreti- 
cally and experimentally, the combined processes of opti- 
cal rectification and the linear electro-optic effect lead to 
an additional, nonresonant contribution into an effective 
nonlinear refractive index of noncentrosymmetric mate- 
rials due to cascading processes. 

The effect of optical rectification is usually neglected 
in the theory of quadratic solitons because the equation 
for the dc field can be integrated explicitly, leading to a 
nonresonant contribution into the effective cubic nonlin- 
earity of the nonlinear Schrodinger (NLS) equation de- 
rived by means of the asymptotic technique in the ap- 
proximation of cascaded nonlinearities (see, e.g., Ref. ff). 
However, for the propagation of spatio-temporal multi- 
dimensional optical pulses in nonresonant quadratic me- 
dia, such a reduction is no longer possible and, as a re- 
sult, the multi-dimensional NLS equation becomes cou- 
pled to a dc field R], similarly to the integrable case of 
the Dawey-Stewartson equation M . 

In this paper we show that the physical situa- 
tion is qualitatively different for periodic quadratically 
nonlinear optical media. We demonstrate that cou- 
pling between the forward and backward waves in one- 
dimensional shallow Bragg gratings with a quadratic non- 
linearity is accompanied by a coupling to the induced dc 
field that appears within the same approximation and 



cannot be eliminated by integration. This effect has been 
overlooked previously, but it leads to a novel physical 
model for gap solitons in quadratic media which we in- 
troduce and analyze here. 

We consider propagation of an optical pulse in a pe- 
riodic medium with a quadratic x^ nonlinear response. 
To derive the coupled-mode equations for the wave enve- 
lope, we start from Maxwell's equation, 
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where V 2 stands for the Laplacian, c is the speed of 
light in vacuum, E is the x-element of the electric field, 
E = E(z,t)e x , and the quadratic nonlinearty is repre- 
sented by a tensor element y' 2 ) = v*- 2 -* • We assume that 
e(z,uj) is a periodic function of z, so it can be presented 
in a general form as a Fourier series, 



e(z,u) = e(u) | 1 + 
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where d = n/k is the period of the Bragg-grating struc- 
ture. Deriving the couple-mode equations below, we as- 
sume the case of a shallow grating, i.e. that the condition 
ej <C 1 holds. Additionally, we may consider a periodic 
modulation of the nonlinear quadratic susceptibility tak- 
ing x^ (z) = x^ ( z + d) as a periodic function with the 
same period d. However, we have verified that this effect 
does not modify qualitatively the analysis and results 
presented below, so that we consider the simplest case 
when x^ 2 -* is constant. 

For a periodic structure, the Bragg reflection leads to 
a strong interaction between the forward and backward 
waves at the Bragg wavenumber ks ~ k. To derive the 
coupled-mode equations for the wave envelopes, we con- 
sider the asymptotic expansion for the electric field in the 
form, 
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where E± = E±{z,t) are slowly varying envelopes of 
the forward(-l-) and bakward(-) waves. The frequency u) 



1 



satisfies the dispersion relation for linear waves, c 2 k 2 = 
Lu 2 e(uj). Due to quadratic nonlinearity, the expansion (|^) 
includes higher-order terms at the frequency 2lj and the 
zero- frequency term, so that the slowly varying functions 
jrj(n,m) _ £j(n,m)^ are defined as nonlinear amplitudes 
of the (n, m)-order harmonics e - mujt e mkz . 

Introducing a small parameter s, we assume E± ~ 
0(e), <9£±/di - aE±/<9z ~ 0(e^ and ej ~ 0(e 2 ). 
Then, substituting the expansion (||), (||) into Eq. (jl]), 
we compare the terms of the same order in front of the 



coefficients e 
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. At the orders (2,0), (2, ±2) and 



(0, ±2) we respectively obtain, 



(2,0) _ 



E 



E (2,±2) = 



2X 



(2) 



-E, E_ 



e(2u>) 



W 2 X (2) 



0( £ 2 ), 



c 2 fc 2 _ w 2 £ ( 2w )] 



0(e 2 ), 



where we have assumed non-resonant interaction with 
the second harmonic, i.e. ui 2 e(2uj) ^ c 2 k 2 . 

At the orders (1,±1) and (0,0) we obtained a system 
of coupled nonlinear equations, 



A\E+\ 2 + B\E„\ 2 + CE(°'°A E+ =0 
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B\E+\ 2 +A\E-\ 2 + CE(°^)E-=0, (5) 
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where v g (uj) = dtv/dk, vq = v g (0), k — w 2 e(w)ei/ (to), 
A = 2( X (2) ) 2 w 4 {/(co)[c 2 fc 2 - uj 2 c(2uj)}}- 1 , B = 
-4( X (2) ) 2 w 2 [/(w)e(2w)]- 1 , C = 2cj 2 x (2) .T 1 H, and 
D = -2x (2) /c 2 with f(ui) = [ui 2 e(ui)y. If we keep the 
transverse coordinates (x,y), Eq. (|6|) should also include 
the transverse Laplacian in the same order. System (|J)- 
(||) describes the interaction between the forward and 
backward waves coupled to a dc wave induced via the 
rectification effect. Importantly, the order of the dc wave 
^(o,o) j g £ 2^ g0 foa.t the dc field itself is of a higher order 
in comparison with the forward and backward scattering 
waves. However, the dc field is coupled to the fields isl- 
and E- in the main order. Importantly, if we assume 
much stronger dc field, e.g. of order of O(l), the system 
is decoupled and the dc wave satisfies an independent 
equation. Including the optical Kerr effect, we obtain 
the same equations as Eqs. ([I]), (|J) with the modified 
constants A and B. 



In the case of a single wave propagating in a homoge- 
neous medium, the induced dc field is explicitly given by 
the host wave || . The similar result is valid for the case of 
a deep grating described by the modulations of the Bloch 
waves, but not for a shallow grating we discuss here. If 
we assume the scaling dE±/dz ~ 0(e 2 ) (as usually done 
in the analysis of higher-dimensional systems such as the 
Dawey-Stewartson equation) and E^ 0,0 ) of order e 4 , then 
the coupling between the dc wave and host-wave scan 
be neglected. However, for isotropic scaling as presented 
here, the effect of the dc wave is directly included into 
Eqs. (||) and (|J). Interaction between the dc field and 
fundamental harmonics has been also discussed in Ref. 
Q, however in that analysis, the dc field appears as a 
cascading effect and its velocity is almost the same as 
the phase velocity. We notice that in our case, the dc 
field is essentially excited by quadratic nonlinearity and 
no assumption is required for the velocity. 

We are looking for spatially localized solutions of Eqs. 
(q)"@ f° r bright gap solitons in the form, 
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where ( = z — Vt; the functions /(£) and #i,2(C)j an d 
the parameters f2, V, A are assumed to be real. The 
parameter g is the argument of the coupling parameter 
k, i.e. k — \n\e lg . Substituting the ansatz ([[) into Eq. 
(pf), we obtain 
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and therefore the contribution of the dc field should van- 
ish at V = 0. 

From Eqs. (|i|) and (||), we set the parameter A as 
yj (v g — V)/(v g + V), and then obtain a system of cou- 
pled equations for /, = £q — 9 2 and 9 + = 9\ + 9 2 , 
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where fi = \k\/(v 2 - V 2 ) 1 / 2 , v = -(2v g n)/(v 2 - V 2 ), 
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We notice that in the problem under consideration, we 
should assume \V\ < v g and \v\ < 2/i. 

Similar to the analysis presented in Ref . |£j] , from Eqs. 
(]§|)-(]l0|) we obtain a closed differential equation for the 
function in the form of the double sine-Gordon (DSG) 
equation. The DSG equation can be integrated, and its 
localized solutions are two types of kinks and anti-kinks 
|Q0|. Using the relevant solutions, we can then find 



/(C) = 
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where the signs ± stand for the cases 5 > and S < 0, 
respectively. Functions 9\ and 82 are then obtained as 



= -2 tan 1 I J— — - tanh T 
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The effective renormalisation of the coefficients due to 
the induced dc field seems extremely important for the 
soliton stability. Indeed, when vq < v g the coefficients 
have a singularity provided V — > vq, changing the char- 
acter of the dependence of the soliton parameters and 
the system conserved quantities on V. The recent stabil- 
ity analysis of the conventional gap solitons |l3| revealed 
the existence of two types of instabilities, oscillatory and 
translational. The most important, translational insta- 
bility appears for large V, so that the induced dc field 
is expected to have a strong effect on the soliton sta- 
bility. In fact, we anticipate that all gap solitons for 
vq < V < v g may become unstable, similar to the soli- 
tary waves in hydrogen-bonded molecular systems [ p^[ . 
It is also worth to mention that in the limit V — *■ vo when 
the coefficients grow, the transverse effect in Eq. (|^) be- 
come important and should be included to compensate 
the singularity. 
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To demonstrate that the issue of the soliton stability 
becomes nontrivial for this model, we present the sys- 
tem invariants. Similar to some other models, we are not 
able to present Eqs. ([|)-@ in a Hamiltonian form di- 
rectly, and therefore we introduce an auxiliary function 
(j) through the relation, ad<j>/dz = E^ ^ - v$D(\E+\ + 
+ C±, (12) l-E-l 2 ), where a 2 = VqD/C. Then, we define the second 
canonical variable as ijj = v^dip/dt and show that Eqs. 
@"@ can b e wl "itt en as a Hamiltonian system 



where C± are integration constants. In order to obtain 
solutions for gap solitons, we restrict the possible an- 
gle variable by the domain, < < 2tt. The solu- 
tion obtained from Eq. ( |l l| ) describes a two-parameter 
family of gap solitons, spatially localized waves in the 
Bragg gratings, which are similar to the gap solitons of 
the conventional couple- mode theory. Actually, by renor- 
malizing the variables as A — > ±tr, B — > ±1, C — > 0, 

1,2 cosQ, 
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C± -> (1/2 ± 1/2)tt ± (4cn>a 2 )/(l - v 2 )ir + 2(f>, we can 
demonstrate that the solution is essentially the same as 
that earlier obtained in Ref. However, the effect of 
the dc wave is included in the parameters 6 and 77. 

In the case \v\ > 2/i, spatially localized solution of Eqs. 
®; @ do not exist. Instead, the kinks of the DSG equa- 
tion for (9_ give solutions for dark gap solitons (see also 
}L2|), localized waves on nonvanishing backgrounds, 
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with the following Hamiltonian, 
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where A = A+v 2 CD and B = B+v 2 CD. Other integrals 
of motion of the system (Q)-@ are the field momentum, 



P 



+ OC 



dz I IE. 



dE; 
dz 



E_- 



dE* 

dz 



dz 



where 




Upper and lower signs correspond to two types of such 
solitons, with the maximum intensity large or smaller 
than the background intensity. 



the total number of the forward and backward waves, 
and an independently conserved number of the dc waves, 
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dz(\E^ 
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dzip. 



Therefore, in a sharp contrast to the conventional theory 
of gap solitons, the model (@)-(q) possesses one additional 
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integral of motion, it has no analogy with other soliton- 
bearing nonintegrable models where the soliton stability 
has been investigated so far. 

In conclusion, we have demonstrated that in peri- 
odic optical media with a quadratic nonlinear response, 
Bragg-grating-induced coupling between the forward and 
backward propagating modes leads always to an induced 
dc field which plays an important role for nonlinear pulse 
propagation in periodic x' 2 ' media. We have derived a 
novel model of the coupled-mode theory for optical gap 
solitons in quadratically nonlinear Bragg gratings, and 
have found the analytical solutions for moving bright and 
dark gap solitons. Analysis of the soliton stability for the 
model (|4|)-(^|) will be presented elsewhere. 
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